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Abstract. We investigate the damping of the baryon acoustic oscillations (BAO) 
signature in the matter power spectrum due to the quasi-nonlinear clustering of density 
perturbations. On the basis of the third order perturbation theory, we construct 
a fitting formula of the damping in an analytic way. This demonstrates that the 
damping is closely related with the growth factor and the amplitude of the matter 
power spectrum. Then, we investigate the feasibility of constraining the growth factor 
through a measurement of the damping of the BAO signature. An extension of our 
formula including higher order corrections of density perturbations is also discussed. 
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1. Introduction 

The baryon acoustic oscillations (BAO) imprinted in the galaxy power spectrum have 
recently attracted remarkable attention, as a useful probe for exploring the origin of 
dark energy [H |2]. The BAO signature has been clearly detected in the 2dFGRS and 
the SDSS galaxy samples P HJ El El [7j. The feasibility of constraining the equation 
of state parameter of dark energy w is demonstrated [H El [10], where w is defined by 
w = Pde/ Pde, where pde and p^e are the pressure and the energy density of the dark 
energy component, respectively. Furthermore, a lot of BAO survey projects are under 
progress, or planned [m [121 [131 [H [ISl [16] . 

Though these BAO surveys will precisely measure the galaxy power spectra, but 
we also need precise theoretical templates, in order to obtain a useful cosmological 
constraint from observational data. In practice, we need to solve uncertainties about 
the gravitational nonlinear clustering of the density perturbations, the redshift-space 
distortions, and the galaxy clustering bias, because these uncertainties might yield some 
systematic effects in the BAO signature. Then, a lot of works related to these topics 
have been done recently, which are very important and challenging for future research 
of dark energy [I3[T8l[l9l[20lini[22l[23l[2l[25l[26l|27^ 

In the present paper, we focus on the damping of the BAO signature in the matter 
power spectrum due to the nonlinear gravitational clustering. We investigate this 
damping in a semi-analytic way on the basis of the third order perturbation theory 
of density fluctuations. The nonlinear clustering is a consequence of mode-couplings of 
the density fluctuations and the peculiar velocity divergence in Fourier space. Especially, 
we demonstrate how the damping of the BAO signature can be expressed with/without 
P22{k) and Pi^{k), which describe the mode-coupling (see sections 2 and 3). Then, we 
develop a simple fitting formula relevant to the damping, which demonstrates the fact 
that the damping of the BAO signature is closely related with the growth factor and 
the amplitude of the matter power spectrum. This suggests that a measurement of the 
damping of the BAO signature might be a probe of the growth factor of the density 
fluctuations, though we need to further investigate if the effects of the redshift-space 
distortion or the galaxy clustering bias are influential to the damping. 

This paper is organized as follows; In section 2, we briefly review the third order 
perturbation theory, and derive the second order power spectrum. In this formalism, the 
nonlinear gravitational clustering is described by the coupling of the Fourier modes of the 
density fluctuations and the peculiar velocity divergences. In section 3, we investigate 
the damping of the BAO signature due to the coupling of the Fourier modes in an 
analytic way. And a fitting formula of the damping, applicable in weakly nonlinear 
regime, is developed. The formula is compared with a result of iV-body simulation. 
Extensions of the formula are also discussed. In section 4, with the use of the fitting 
formula, we demonstrate a future feasibility of constraining the growth factor through 
a measurement of the damping of the BAO signature. The last section is devoted to a 
summary and conclusions. Throughout this paper, we use the unit in which the velocity 
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of light equals 1, and adopt the Hubble parameter Hq = 100/ikm/s/Mpc with h = 0.7, 
unless explicitly stated. 

2. Quasi-nonlinear evolution of the matter power spectrum 

To investigate the effect of the nonlinear gravitational clustering in the present work, 
we employ the standard perturbation theory (SPT) of the density fluctuations and the 
peculiar velocity divergences [281 ISH EOl EH |32l [33l |Ml ESI ES] . Numerical simulation 
is a rigorous approach to the gravitational clustering, but the perturbative approach 
has an advantage to understand which factors are relevant to the damping of the BAO 
signature in an analytic way. Recently, several authors have developed new formalism 
beyond the SPT ^EHlE9lllQlliTlll2lll3lllllll5lll6]. In the later subsection 3.3, we 
discuss a possible extension of our main result to include these formalism. However, let 
us start with briefly reviewing the derivation of the second order power spectrum on the 
basis of the SPT up to the third order of perturbations. 

We consider the matter fluctuations after the recombination whose wavelength of 
interest is smaller than the horizon size, then the evolution of the matter fluctuations 
can be analyzed by the pressure-less nonrelativistic fluid with the Newtonian gravity. 
Denoting the comoving coordinates by x, and the conformal time by rj, the evolution 
equations are 

5ix, r^) + V ■ [(1 + 6ix, 7])) v{x, t])] = 0, (1) 
v{x, 1]) + {v{x, 7]) ■ V) v{x, 7]) + H{r])v{x, r]) = -V(j){x, t]), (2) 

v'<p{x,v) = ^n\v)H^,v), (3) 

where 6 is the density contrast, v is the peculiar velocity, is the gravitational potential, 
a is the scale factor, a dot denotes the derivative with respect to the conformal time, 
and Ti = a/a. Here, the Einstein-de Sitter universe is assumed. 

We ignore the rotational mode of the velocity, since our interest is only the growing 
solution, and the rotational mode is the decaying solution in the expanding universe. 
Then, we introduce the velocity divergence, 

9{x,r]) = V ■ v{x,r]). (4) 

It is convenient to analyze the perturbations in the Fourier space, and we define the 
Fourier coefficients as 

Six,r^)= [ -^5(fc,ry)e-'»-, (5) 



0{x,v)= I -—eik,v)e-"'''. (6) 



(27r)3 

Then, equations ([I]) - (jS]) yield 



S{k, 7]) + e{k, V) = - j d^'ki I (fk25^^\ki + k2- k)^^e{ki,T])6ik2, r]), (7) 



(8) 



where 6^^\k) denotes the Dirac's delta function. The right hand side of these equations 
describes the mode-couphngs which govern the nonhnear evolution of the matter 
fluctuations. 

To solve these coupled equations, we adopt the perturbative expansion as 



S{k,7^) = ^a"(r;)5,(fc), e{k,7j) = 7i(r;) ^a"(r/)^n(fe). 

n=l n=l 

In general, the n-th order solution can be written as 

Sn{k) = J d\ • • • y d'qj^'^ - ^) ^n(9l, ...,qn) R ^l(li)^ 

Onik) = -J d% ■ • • I d^qj^^^ (^XJ - ^ Gniqu ...,qn) fl S^, 
where . . . , q„) and Gn{qi, ■ ■ ■ , qn) are defined as 



(9) 

(10) 
(11) 



i;, / \ 'y^^ . . . ,q 

Fn{qi,...,qn) = 2^ 

m 



\ (2n + 3)(n- 1) 

k\k^ ■ k2) 



k ■ ki 

1 + 2n) Fn-m{qm+l, ■ ■ ■ , qn) 

Ki 



+ 



n i ^ ^ ^ _ \ " ^miqi, ■ ■ ■ 1 qm) 

^n[qi: ■ ■ ■ iqn) — / 



k ■ ki 



\ (2n + 3)(n- 1) f k\ 
k'^{ki ■ k2] 



G n—m{,Qm+l 1 ■ ■ ■ t qn) 



-Fn—m{,qm+lj • • • ) 9n) 



(12) 



+n- 



2u2 



'Gn—m{^Q.ra+li • • • i Qn) 



(13) 



Assuming that the first order density perturbations described by 5i{k) is a Gaussian 
random field, we obtain the second order matter power spectrum 

^spt(A;, z) = Dl{z)PUk) + Dl{z)P2{k), (14) 
where P\in{k) is the linear power spectrum given by 

(27r)35(3)(fc + k')PUk) = {5,ik)5,ik')), (15) 

Di{z) is the linear growth factor, and P2{k) is the second order contribution to the 
power spectrum, which is conventionally expressed as follows; 

P2{k)^P22ik) + 2P^sik). (16) 

Taking the 4-point correlations of Si{k) into consideration, we obtain the explicit form 
of P22ik) expressed as the integral of the square of the linear power spectrum, 

P22{k) = 2 J d'qPi^{q)Pi^ {\k-q\) [F|(q, k - q)f , (17) 



5 



where ^2) is symmetrized over its arguments. On the other hand, Pi3{k) has the 

form shghtly different with ^22(^)5 

2P,,{k) = 6P,Uk) J d^qPnn{q)F^{q, -q, k). (18) 

The sohd curve in figure 1 shows typical behaviour of P22{k) and 2Pi3(A;), where the 
cosmological parameters are h = 0.7, flm = 0.28, Qb = 0.046, Ug = 0.96 and as = 0.82. 
We discuss further details of the typical behaviour of P22{k) and 2Pi^{k) in the next 
section. 

Originally, this perturbation expansion is consistently formulated in the Einstein 
de-Sitter universe. The result is extensively used for the other general cosmological 
model, with the corresponding growth factor normalized as Di{z) = a{z) at a{z) <^ 1. 
We follow this prescription because the validity of the extensive use of the result is 
known, for example, in the literature [361 W7\ . 



3. Damping of the BAG signature 

In this section, we will examine the damping of the BAO signature due to the nonlinear 
gravitational clustering with employing the third order perturbation theory. The BAO 
signature in the matter power spectrum can be extracted as follows: 

(19) 

where P{k, z) is the matter power spectrum including the BAO signature, but P{k, z) 
is the matter power spectrum without the BAO, which is calculated using the no-wiggle 
transfer function in [18]. Hereafter, the quantity with the 'tilde' implies the quantity 
computed using the no-wiggle transfer function. Within the third order perturbation 
theory, we may write 

PsPT(fc, z) = Dl{z)PUk) + D{{z)P2{k), (20) 
where P2{k) is 

P2(A:) =P22(fc) + 2Pi3(fc), (21) 

and P22{k) and Pi^ik) are, respectively, defined by ( |T71) and ( |T8l) . but with the no-wiggle 
transfer function. 

Figure 2 shows B^^^{k, z) as a function of the wavenumber k for several redshifts, 
which is obtained using the third order perturbation theory (see also below). Here, the 
spatially flat universe with the cold dark matter (CDM) and the cosmological constant 
A is assumed, where the cosmological parameters are the same as those of Figure 1. 
The oscillating behaviour of the curves is the BAO signature. As the redshift becomes 
small, one can see that the amplitude of the oscillations decreases. This damping of the 
oscillations is more significant as k is larger. 
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3.1. Analytic approach 

The aim of this paper is to understand the nature of the damping of the BAO signature 
in detail. To this end, the analytic approach based on the third order perturbation 
theory is useful. With the use of the formula of the second order power spectra, 
and ( |20l) . we have 

z) = ^^P^^^'^^ - 1 = PUk) + Dl{z)P2{k) _ ^ 
'^^^'^ Pspt(A;,^) PUk) + Dl{z)P,{k) ■ ^ ' 

First, we adopt an approximation, 

P2(fc)~P22(fc)+2Pi3(A;). (23) 

In brief, P22{k) in f|T6l) is replaced with P22{k)- The validity of this approximation is 
demonstrated in figure 1, where the upper solid curve is ^22(^)5 and the upper dotted 
curve is P22{k)- The validity of this approximation comes from the fact that the mode- 
coupling of different Fourier modes decrease the coherent BAO signature. As shown in 
this figure, the tiny oscillatory feature remains. This may induce somewhat the small 
shift of the peaks(troughs) of the BAO, as mentioned in the reference [H]. However, 
this shift of peak location would not be problematic, as long as we are interested in the 
damping of the BAO signature. 

On the other hand, Pi3{k) can not be simply replaced with Pi^lk). However, careful 
consideration leads to an expression for Pi3{k) in terms of the linear power spectrum 
multiplied by a monotonically decreasing function, as follows. First, we define 

BUt) . - 1. (24) 

which corresponds to fl22|) . within the linear theory of density fluctuations. Note that 
i?iin(/c) is not time-dependent. With this definition, we obtain 

2Pi3(A;) = 6Ph,(A;) [1 + B^k)] J d^qPiMF^i^, "Q, k). (25) 

Here, we apply the following approximation to the linear power spectrum of the 
integrand, 

2Pi3(A;) ~ QPy^k) [1 + BUk)] J rf^gPun(g)i^3 (q, -q, k) 

= 2[l+Pu„(A;)]Pi3(A;). (26) 
Substituting fl2Bl) into (1251) . P2{k) is written as 

P2{k) = P2{k) + 2BUk)Piz{k). (27) 
Then, from (1221) . we obtain 

1 + D?(.)'^-^') 



PrPT*(fc, z) ^^BUk). (28) 

^PiUk) 
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This formula indicates how the BAO signature is modified as the gravitational clustering 
evolves, which is expressed by the BAO signature in the linear theory multiplied by the 
correction determined by the no-wiggle quantities and the growth factor. The second 
term of the denominator in fl28|) is small within the perturbation scheme, we may expand 
it as 




PsPT {k, z) 

, -..v ., (29) 
PUk)) J J 

though the higher order terms make no sense because we are working in the second 
order theory of the power spectrum. The expression (129|) indicates that the leading 
effect of the nonlinear mode-coupling on the damping is described by the factor 
-Dl{z)P22{k)/Piin{k), and that the sign of the term clearly shows that this effect is 
a damping. 

3.2. Fitting formula 

As we are working within the third order perturbation theory, its prediction does not 
perfectly coincide with the result of full order computation, which can be obtained 
by A^-body simulations. However, we believe that the prediction of the third order 
perturbation theory is useful in constructing a semi-analytic formula which reproduces 
the result of A^-body simulations. Then, we here find a simple fitting formula which 
reproduces the prediction of the third order perturbation theory. We discuss the validity 
of the formula in comparison with a A^-body simulation in the below. 
Up to the second order of Di{z), fl29l) yields 



-^SPT {k, z) 



.2,.^^22(A:) 



BUk)- (30) 




(31) 



Piin(A;). 

From a detailed analysis of P22{k) / P\m{k) as a function of /c, we find that the following 
fitting formula works well, 

P22{k) " 

PrUk) 

where ag is the rms matter density fiuctuations averaged over the sphere with the radius 
of 8/i~^Mpc, kn and 7 are the constant parameters which depend on flmh"^, ^bh^ and 
Hg. Figures 3 and 4 show the best-fit value of kn and 7 as a function of VLmh'^ and 
where we fixed Hs = 0.96. This shows that these two parameters depend on Qmh"^ and 
flbh'^ linearly. We can show a similar dependence on Ug. Then, we found the following 
fitting formula: 

fc„ = -1.03(fi„/i^ + 0.077)(r];,/i^-0.24)(n, + 0.92) /iMpc"\ (32) 
7 = -11.4(fi„/i2- 0.050)(^]^,/^^- 0.076) (n,- 0.34) hMpc'\ (33) 
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Though these dependence on the cosmological parameter might have to be investigated 
more carefully, but the validity is guaranteed in the following narrow range 

0.13 < rimh^ < 0.15, (34) 
0.022 < ribh^ < 0.024, (35) 
0.94 <n,< 0.98. (36) 

In the present paper, we simply assume that the dependence on the other cosmological 
parameters except for flmh"^, ^bh^ and Ug can be ignored. 

Finally, we obtain the heuristic expression of the leading correction for the BAO 
signature. 



Bf^^{k,z)= I- f{n^K',n,K\ns,a^,k,z) BUk), (37) 



with 



f{n^h:\nbh:\ns,a8,k,z) = aiDt{z){-^\ (1-7). (38) 




k 



Figure 5 demonstrates the agreement between this fitting formula, -BgpT, and the 
prediction of the third order perturbation theory, -Bgpx*, (l22l) . As one can see from 
this figure, the agreement becomes worse as the redshift is lower and the wavenumber is 
larger. Figure 6 shows the relative error of the fitting formula, - 5§pf *|/|5§pt*|, 

as a function of the redshift at the wavenumbers of PI, P2, P3, Tl, T2 and T3, which 
are defined in Figure 5. The relative error is about less than 10 % at worst in the range 
of the wavenumber k ^ 0.2 /iMpc~^ until z ~ 1. The sign of -Bfpx ~ -Bgpx* of the third 
peak (trough) changes around the redshift 2.7 (1.4), where the relative error becomes 
zero. For the redshift less than 1, the agreement becomes worse especially at the large 
wavenumber k ( ^ 0.2 hMpc~^). However, we should also note that the third order 
perturbation theory becomes worse to reproduce iV-body simulations for lower redshift 
and for larger wavenumbers. 

Figure 7 shows a comparison of our fitting formula {Bg^prp, solid curve) and a result 
of A^-body simulation (squares with error bar) [19]. This demonstrates that our fitting 
formula reproduces the result of A^-body simulation for k ^ 0.2 /iMpc~^ until 2; ~ 1 
within error bars, roughly. In this conclusion, note that we focus on the damping of the 
BAO signature, not on the amplitude of the power spectrum itself. 

In the A^-body simulation, we adopt a ACDM cosmology {Qm = 0.279, fl\ = 0.721, 
nh/n„, = 0.165, h = 0.701, Us = 0.96, ag = 0.817; WMAP5 best fit value, [50]), and 
calculate the linear matter power spectrum using CAMB ^51j. We adopt 512^ particles 
in periodic cubes with each side 1000/i^^Mpc, and displace A^-body particles using the 
second-order Lagrangian perturbation theory (e.g., [52]) from uniform grid positions at 
z = 31. The simulations are carried out using the Gadget2 code [53] to output data of 4 
redshifts [z = 3, 2, 1, and 0.5). Our total simulation volume is 8/i~^Gpc^, which might 
be small to investigate our scales of interest (first a few BAO peaks). Then we correct 
the deviations from the ideal case of infinite volume, as follows. 
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In addition to the time integration using Gadget2, we also calculate the time 
evolution of the density contrast using the second-order perturbation theory starting 
from the identical initial condition with that used for the A^-body simulation. We then 
calculate the power spectrum from the result of the perturbation theory, and measure 
the deviation from the the linear power spectrum at each output time. The deviation 
of our finite volume simulation is well explained by the mode-couplings predicted by 
the second-order perturbation theory. Then, we can obtain a corrected power spectrum 
by multiplying the measured power spectrum from A^-body simulations by the ratio of 
the linear power spectrum to the power spectrum predicted by the perturbation theory. 
Thus the error bars in the figure stand for the remaining standard errors after this 
correction, which are very small compared with the usual cases in the region of the 
small wavenumbers. This is because the perturbation theory used for the correction is 
more accurate at the small wavenumbers, which can remove the deviations from the 
linear power spectrum better. This correction does not improve the errors at large 
wavenumbers (see |2l], [19] for more details). Moreover, in order to construct the no- 
wiggle spectrum required for computing the ratio B, we use a cubic basis-spline fitting 
with the same break points as in [27] (see also [8]). 

3.3. Possible extension and discussion 

There are other possible effects which may affect the damping of the BAO signature: 
the higher order nonlinear effect, the redshift-space distortions and the clustering bias. 
These effects may be influential to the damping of the BAO signature, however, there 
remains a lot of uncertainties about these effects, at present. Here, let us consider a 
possible extension of our work to include the higher order nonlinear corrections in real 
space. 

The standard perturbation theory is useful to analyse the damping of the BAO 
signature in an analytic way, however, it is not enough for precise predictions that match 
with result of A^-body simulations in the regime where the nonlinear effect becomes 
significant. Recently, several authors have developed non-perturbative approach to the 
nonlinear density clustering, as mentioned above [87l IM IM 1401 14T1 l42l l48l l44l l45l l46] . 

As an alternative to the SPT, we here consider the work proposed by [36], which 
uses the technique of resumming infinite series of higher order perturbations on the basis 
of the Lagrangian perturbation theory (LPT). One of the advantage of this approach is 
the simplicity of the resulting expression of the nonlinear power spectrum, which enable 
us to incorporate the result into our formula. In the framework of the LPT|l6], the 
matter power spectrum can be given by 



PhPT^k, z) = e 



(^)'^(^)[Di(^)2pHn(fc) 



(39) 



where 




(40) 
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Corresponding to 
based on the LPT. 



we define the BAO signature of the matter power spectrum 



LPT 



{k,z) 



- 1, 



(41) 



PLPT{k,z) 

where Plpt(^5 z) is defined by fl39|) but with the no-wiggle transfer function. Adopting 



an approximation, g{k) 



BrFT\k,z) 



g{k), we can obtain the following expression 
Dliz) 2Pi3(A:) 

hn{k). 



1 + Dj{zyg{k) Ph„(A:) 



Djiz) 



P2{k) 



By, 



(42) 



1 + Dj{zrg{k) PiUk) 

Repeating the same procedure from (125]) to (|5U]) . one can obtain the leading correction 
to the BAO 



B^PTik^z) 



Dl{z)~g{k) P^^{k) 



BUk), 



then we obtain 



Slpt(^; z) 



1 



f{nmh'^,nbh'^,ns,as,k,z) 



BUk), 



(43) 



(44) 



1 + Df{zrgik) 

as an extended fitting formula. Note that this reduces to the expression fl3UI) in the limit 
of g{k) — >• 0. Comparing (jH]) and fl37j) . the difference is the contribution from Df{a)g{k) 
in the denominator in front of P22{k) / P\in{k). Since g{k) is positive, this correction make 
the damping of the BAO signature weaker compared with (1371) . In figure 7, the dashed 
curve plots Bl^rj,{k, z), (l44l) . 

In order to show the validity of the fitting formula, figure 8 plots the relative error 
s \6B/B\ = I-Blpt ~ -^lFt* I / 1 -^LPT* I wavenumber of the peaks and troughs as a 

function of the redshift, which is the same as Figure 6 but with the LPT instead of the 
SPT. As one can see from this figure, the fitting function -Blpt better reproduces the 
damping of the BAO signature -Blpt*; compared with the case of the SPT, especially 
for larger wavenumber even at lower redshift. The relative error is less than 10 % in the 
range of the wavenumber k ^ 0.2 /iMpc~^ (near the third peak) until the present epoch, 
z = 0. In the below, we investigate other possible ways in describing the damping of 
the BAO signature due to the nonlinear effect. 

One of the other possible formulas for the fitting function is the exponential 
function, as has been discussed in references (e.g., [71 [20l |22l [M])- Angulo et al [22] 
discussed the following fitting formula with the Gaussian damping function, 

, 21 



B{k^ z) = exp 



k 



.V^k,{z), 



Bn^{k,z) 



(45) 



where fc*(z) is a time- dependent free parameter, which should be calibrated by 
measurements from A^-body simulations. This free parameter describes the time- 
dependence of the BAO damping, and Sanchez et al [51] obtained the best fit value. 
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k^{z) = 0.172 /iMpc~^ at z = 1, from their result of A^-body simulations in Their 
model includes another parameter, a, which describes the shift of the BAO scale in the 
wavenumber. However, in the present paper, we adopt the case of no shift of the BAO 
scale (a = 1) because our interest is focused on the damping of the BAO signature not 
on the shift of the BAO scale. It is also shown that a almost equals 1 [21] . 

Figure 9 compares our fitting functions and the Gaussian damping function. The 
vertical axis stands for B^^{k, z = 1) divided by B\[^{k, z = 1), which means the damping 
function. The dotted curve and the sohd curve are -BfpT and -Blpt^ respectively. The 
dashed curve is the Gaussian damping function, (HSl) . with k^{z) = 0.172 /iMpc~^. -BfpT 
becomes negative as the wavenumber becomes larger, k ^ 0.25 /iMpc~^. This is because 
the phase inversion of the BAO signature appears in the SPT, as shown in the panel of 
z = 1 in Figure 5, where of the SPT breaks down in this regime as mentioned in Section 
3.2. On the other hand, the damping function based on the LPT, i?Lprp/5iin, shows the 
similar behavior as the Gaussian damping function, which approaches zero in the limit 
of — > 0.7. This is due to the modification by the factor, 1/[1 + Df{z)g{k)], in (jSj). 
-Blpt/Aih becomes negative for k ^ 0.7 hMpc^^ because of the same reason as that for 

In figure 9, we also plot the following damping function (dot-dashed curve), 

Bll\{k, z) = exp [-/(a^/i^ ribh'^, n^, as, k, z)\ Bi^^{k). (46) 

The leading term of the expansion of this damping function with respect to D\{z) leads 
to 5f^T. (inZD- As one can see, this damping function -Bgpx agrees with the Gaussian 
damping function, ( l45l) . This means that the higher order perturbations is important 
in describing the BAO signature for the regime of the wavenumber, k ^ 0.2/iMpc~^ at 
redshift 1. 

The result of this section gives us a clue to find how to describe the BAO damping 
due to the nonlinear gravitational clustering. The BAO damping is determined by the 
normalization as and the growth factor Di{z). This suggests that a measurement of the 
BAO damping might be useful to estimate the growth factor of the amplitude of the 
density perturbations, asDi{z). In section 4, we demonstrate a feasibility of constraining 
asDi{z) by measuring the damping of the BAO signature. 



4. Feasibility of Constraining asDi{z) 

As shown in the previous section, the damping of the BAO signature is closely related 
with the amplitude of the power spectrum, which is determined by asDi{z), while the 
BAO signature within the linear theory is determined by the density parameters flmh"^ 
and Qfih"^. In this section, we discuss the feasibility of constraining asDi{z) by measuring 
the damping of the BAO in the power spectrum in quasi-nonlinear regime. To this end, 
the formula developed in the previous section is useful. 

As mentioned in section 1 and subsection 3.3, the redshift-space distortions and the 
clustering bias might be additionally influential to the damping of the BAO signature. 
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However, we here assume an optimistic case that the damping is determined by the 
quasi-nonhnear clustering effect and neglect the effects on the damping of the BAO 
signature from the redshift-space distortions and the clustering bias. Then, we study 
how a measurement of the damping is useful to determine asDi^z). Very recently, it is 
recognized that a measurement of the growth factor of the density perturbations is a 
key to distinguish between the dark energy model and modified gravity model for the 
cosmic accelerated expansion (e.g., [551 156], and references therein). Our investigation is 
the first step to investigate if a measurement of the BAO damping is useful to measure 
the growth factor. 

In our investigation, we adopt a simple Monte Carlo simulation of the galaxy power 
spectrum assuming the ACDM model. The error of the galaxy power spectrum depends 
on its amplitude. For definiteness, we assume the galaxy power spectrum Pggi{k,z) is 
modeled with the no- wiggle linear power spectrum in real space, Piin(fc,2;), as 

P,.,,{k, z) = [l + B{k, z)]P,^i{k, z), (47) 

with 

where we use f l44l) as B{k,z), 6 is a constant bias factor, and Ai, and Q are the 
parameters, which describe the correction of the nonlinear clustering, the redshift-space 
distortions and the scale-dependent bias to the no-wiggle power spectrum. This model 
is based on the Q-model of Cole tt al [5^, and is elaborated by Sanchez et al [51] by 
adding the new parameter A2{= Q/IO) which better reproduces the nonlinear power 
spectrum at large wavenumber. Cole et al [57] showed that, from numerical study, the 
value of Al = 1.4 is adequate to reconstruct the galaxy power spectrum in redshift- 
space, though Q strongly depends on the galaxy type [571 EHl EH] . Then, for simplicity, 
we here adopt Ai = 1.4 and consider the cases Q = 4, 8, 16, 32, to estimate the error 
for the galaxy power spectrum. We assume the bias parameter b = 2, for simplicity. 

To estimate the constraint on the growth factor, we perform a simple Monte-Carlo 
simulation. We assume that the BAO signature, B{k, z), can be extracted from a galaxy 
power spectrum by the method like in the reference [5|. The variance of the error in 
measuring the BAO signature can be estimated by 

APgii(fc,z) 
with 

APl,{k,z) = 2^^Q'{k,z), (50) 

where AV^ = Airk'^Ak is the volume of the shell in the Fourier space, and Q'^{k,z) is 
defined as 

Q-'ik,z) = Aa[ dz^^s'lz]- , r (51) 

^ ^ J dz ^ ^ l+n sz])P^.Jk,szW ^ ^ 



AB\k) = ^/""'^"'"^ (49) 
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where n{s) is the comoving mean number density, s = s[z] is the comoving distance- 
redshift relation, and AA is the survey area. 

With the use of the above formulas, we assess defined by 



where B{ki, zY^ is the theoretical one of a fiducial target model at the wavenumber ki, 
while B{ki,z)°^^ is the corresponding observational one. The fiducial target model is 
Qm = 0.28 and ag = 0.82, Qi, = 0.046, h = 0.7, and = 0.96, which are the same 
as those adopted in figure 1. B{ki, z)"'"' is obtained through a Monte Carlo simulation, 
following the steps in the below, 

(i) Based on the fiducial model, compute B{ki,z) and the variance AB{ki,z) at 
ki = Ak{i — 0.5), for i = 4, 5, ■ ■ ■ , 19. Here we specify a bin of the Fourier space, 
Ak = O.Ol/iMpc"^, and consider the range of 0.03 < k < 0.19, where the validity of 
our formula in the previous section is guaranteed. We assume two galaxy redshift 
samples as typical future survey. One is the WFMOS-like sample, AA = 2000 deg^ 
in the range of redshift 0.5 < z < 1.3, and n = 5.0 x 10""* [/i^^Mpc]~^, which 
contains 2.1 x 10^ galaxies. The other sample assumes the same number density 
and the range of the redshift, but the larger survey area A A = An steradian. 

(ii) Each B{ki,zy^'^ is obtained through a random process assuming the Gaussian 
distribution function with the variance AB'^{ki, z). The data points in figure 10 
show an example of a set of B{ki, generated through the random process. 

(iii) We assess the values of with this sample obtained by the step (i) and (ii). 

We iterate these steps and compute 1000 sets of x^, then obtain the average of x^- 
For the theoretical model, B{ki,zy^, we fixed Qbh"^, h, and as those of the fiducial 
model, but took Qmh"^ and a^Di^z = 0.9) as the variable parameters. Solid curves in 
figure 11 show the contour of Ax^ = 2.3 (inner curve) and Ax^ = 6.17 (outer curve), 
which corresponds to the 1 a and the 2 a statistical confidence level, respectively, in the 
Qmh'^ and asDiiz = 0.9) plane, for the WFMOS-like galaxy sample of AA = 2000 degl 
The dotted curves show the same but with the sample of the survey area AA = Air 
steradian. 

From figure 11, one can read that the 1 a error of asDi{z = 0.9) is about 0.2 for 
the sample of AA = 2000 deg^, and is about 0.05 for the sample of AA = An steradian. 
Thus, the constraint on a^Di^z = 0.9) of the WFMOS-like sample {AA = 2000 deg^) is 
not very stringent. We already have the small uncertainty of asDi{z = 0.9) at the level 
of a few X 0.01 [36j, from the result by the WMAP observations, which is obtained on 
the basis of the fiat CDM cosmological model with the cosmological constant. However, 
we note that the BAO damping is unique and independent to obtain a constraint on 
asDi{z) around the redshift 1. 

It would be useful to discuss the origin of the error. To this end, we adopt the 
Fisher matrix approach. As we are considering the constraint from the damping of the 



[B{ki, zY'' - B{ki, z) 



obs 



(52) 



AB^h.z) 
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BAO signature, we may work with the formula for the Fisher matrix 



~ 47r2 de\ de] Q\k,z) ' ^^^^ 

where 9i denotes cosmological parameters, for which we focus on a^Di{z) and Qmh'^, 
Pg^i{k,z) is defined by (gH]), Q'^{k,z) is defined by ([5l]), we adopt (jSD as B{k,z), 
and fcmin = 0.02/iMpc"^ and fc^ax = 0.2/iMpc"\ The minimum error attainable 

— 1/2 

on is expressed by the diagonal part of the inverse Fisher-matrix, A6i = Fg.g. , 
if the other parameters are known. Figure 12 plots the errors A{a^Di{z = 0.9)) 
and A(f2m^^)) obtained from the Fisher-matrix, as a function of the mean number 
density n(6/2)^lil The thick (thin) solid curve is A{asDi{z = 0.9)) of the sample with 
A A = 2000 deg.^(AA = Att radian). The dotted curve is A(f2m^^)- Here we adopted 
Q = 16, but the result is not sensitive to this choice. One can read that the value of the 
error at the point n{b/2)'^ = 5 x lO^^f/i^^Mpc]^'^ is consistent with the result of figure 
11. 

Figure 12 demonstrates that the constraint can be improved by increasing 
the number density of the galaxy sample, but can not be improved for 
n(6/2)^ ^ 10~^[/i~^Mpc]~"^. The reason can be explained using figure 13, which plots 
nPgai{k,z = 0.9) as a function of the wavenumber, where we fixed n{b/2)'^ = 5 x 
10~^[/i"^Mpc]~^. In the region k < 0.2 /iMpc~\ we have nPgai{k,z = 0.9) > 1, which 
means that the shotnoise is not the dominant component of the error. However, the 
constraint on the parameter is slightly improved by increasing the number density fz, 
depending on the bias b. 

Finally in this section, we mention the effect of the redshift-space distortions and the 
clustering bias on the BAO damping, which might be influential in measuring asDi{z). 
It would be true that there still remains room to investigate how the redshift-space 
distortions and the clustering bias affect the BAO damping. If the effect of the redshift- 
space distortions and the clustering bias on the BAO damping could not be clarified, 
it would make it difficult to conclude that the BAO damping is useful. However, very 
recently, several authors have discussed the issue [Ml EOl |6T] . If the effect on the BAO 
damping are well understood, it can be a signal, and would be useful in measuring the 
growth factor. 



5. Summary and Conclusions 

In the present paper, we examined the effect of the nonlinear gravitational clustering 
on the BAO signature in the matter power spectrum. In particular, we focused on the 
damping of the BAO signature in the quasi-nonlinear regime. Our approach is based 
on the third order perturbation theory of the matter fluctuations, which enables us to 
investigate the damping in an analytic way. We found a simple analytic expression that 
describes the damping of the BAO signature, which clarifies what the important factor is 

I Note that the Fisher-matrix depends on the parameter nb^. 
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for the damping. We showed that the leading correction for the damping is in proportion 
to the combination of {crsDi{z))^. On the basis of the result, we constructed a fitting 
formula for the correction of the damping of the BAO signature in the weakly nonlinear 
regime, which is expressed as a function of k, Qmh'^, Qbh"^ and Ug. This fitting formula 
reproduces the damping of the BAO signature of the second order power spectrum 
within the standard perturbation theory at 10% level for k ^ 0.2 hMpc~^ until 2; ~ 1, 
though the formula is only guaranteed in a narrow range of the parameters flmh"^, ^bh"^ 
and Ug. 

We also discussed a possible extension of our formula to elaborate higher order 
nonlinear corrections using a technique of resumming infinite series of higher order 
perturbations on the basis of the Lagrangian perturbation theory. This extended formula 
reproduces the damping of the BAO signature of the second order power spectrum based 
on the Lagrangian perturbation theory at 10% level for k ^ 0.2 hMpc~^ until 2 ~ 0. 
This formula was compared with a result of A^-body simulation, which showed the 
validity of the extended formula. 

A measurement of the damping of the BAO signature might be useful as a probe 
of the growth factor of the density fluctuations. As a first step to investigate such the 
possibility, we assessed the feasibility of constraining asDi{z) by measuring the damping 
of the BAO in the power spectrum. For a useful constraint, we need a very wide survey 
area of the sky. For a definite conclusion, however, we must include the effect of the 
redshift-space distortions and the galaxy clustering bias on the damping of the BAO 
signature. Thus, more sophisticated formula including the effect of the redshift-space 
distortions and the galaxy clustering bias at the same time is required. If this effect on 
the BAO damping was well understood, it can be useful in measuring the growth factor. 
We plan to revisit this issue in future work. 
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Figure 1. Typical beliavior of P22{k) (upper solid curve) and 2Pi3(fc) (lower solid 
curve), which are calculated using the transfer function suggested by [48J with the 
BAO. The cosmological parameters are h — 0.7, flm = 0.28, ilb = 0.046, Us — 0.96 
and as — 0.82. The dotted curves show P22{k) (upper curve) and 2Pi3(fc) (lower 
curve) . 
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Figure 2. i?gp'!f*(fc, z) as a function of k for several redshifts, z — 2, 1, 0.5, which are 
derived from the matter power spectrum including the second order contributions. The 
solid curve is the linear theory. The cosmological parameters are the same as those 
of figure 1. One can see the damping of the amplitude of the BAO as the redshift 
becomes small. In addition, this damping is more significant as the wavenumber k is 
larger. 
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Figure 4. The best-fit 7 as a function of flmh"^ and The other cosmological 

parameters are the same as those of figure 1. 



20 



N 
PQ 



0.1 
0.05 


-0.05 




0.01 
OQ 0.005 



r 



-0.005 
0.1 



0.05 - 



OQ 



-0.05 - 



0.01 
OQ 0.005 



PI 



I— I— p-i— I— I— 

SPT, exact 

SPT. fit 

linear 



t t T3 

Tl T2 
H — I — I — I — I — I — I — I — I — I — I — I — I — h 



- H — I — I — I — — I — I — I — I — — I — I — I — 



H — I — I — I — — I — I — I — I — — I — I — I — h 



SPT. exact 

SPT. fit 

linear 



\ 

-0.005 I-. ■ ■ ■ I 



— I— |— I— I— I— 1 

SPT, exact 

SPT, fit 

linear 



z=2 



H — I — I — I — I — I — I — I — I — I — I — I — I — h 

SPT, exact 

SPT, fit 

linear 




0.1 0.2 

k[hMpc-i] 



0.3 0.1 0.2 

k[hMpc-i] 



0.3 



Figure 5. This figure compares our fitting formula {Bg^prp, daslied curve), (|37p . 
and the formula of the third order perturbation (i?gpx*: solid curve), ((22)l . for various 
redshifts (z=3, 2, 1, 0.5). The dotted curve is the linear theory. 6B is the relative 
difference between the solid curve and the dotted curve {SB — i?gpx ~ -^spt*)- Here, 
the cosmological parameters are the same as those of Figure 1. PI, P2 and P3 (Tl, T2 
and T3) denote the first, the second and the third peak (trough), respectively, used in 
Figure 6. 
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Figure 6. Relative error \5B/B\ = \B^^^ - S^ptN/I-Sspt*! as a function of the 
redshift at the wavenumber of PI, P2, P3, Tl, T2 and T3, respectively, defined in 
Figure 5. 
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Figure 7. The square with the error bar is the result of A^-body simulation. The solid 
curve is the fitting formula based on the SPT, i?gpx of ([37)) . while the dotted curve is 
the linear theory. The dashed curve is the result of an extended fitting formula, Slpt 
of (|44|) . The degree of the damping of -Blpxi is slightly weaker than that of i?gpx- 
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Figure 8. Same as Figure 6 but witli tlie LPT instead of tlie SPT. Relative error is 

\6B/B\ = \Bt\,r^-Blf^'\ 
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Figure 9. Comparison of fitting functions of tlic damping. The dotted curve is i?spT> 
([57)) . while the solid curve is BIpt, (m. The dashed curve is the Gaussian damping 
function, (gSl), with K{z) = 0.172 hMpc~^. The dot-dashed curve is B^^^,, (gS]). In 
this plot, we adopted the same cosmological parameters as those in 
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Figure 10. The solid curve is the prediction of the fiducial model, whose parameters 
are the same as those of figure 1, while the squares with the error bar is an example 

of B°'^"(fcj,z = 0.9), obtained through our Monte Carlo simulation. For the galaxy 
sample, we assumed the WFMOS-like sample of the comoving mean number density 
n = 5.0 X 10~^ [/i~^Mpc]~^ and the survey area 2000 deg^ in the range of redshift 
0.5 < ^; < 1.3. The bin size of the Fourier space is Ak = 0.01 /iMpc~^. 
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Figure 11. Contours of Ax^ in Qmh^ and asDi{z = 0.9) plane. The solid curves 
assume the WFMOS-like sample of the survey area, A A = 2000 deg^. Inner (outer) 
curve is the contour of Ax^ — 2.3 (Ax^ = 6.17), which corresponds to 1 a (2 a) 
confidence level. The dotted curve is the same as the solid curve, but assumes the 
survey area, AA = An steradian. 
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Figure 12. The Icr-level statistical errors of ^m)^ (solid curve) and a^D\{z = 0.9) 
(dotted line) as a function of the number density n{h/2Y. The thick curves assume 
A A = 2000 deg^, and the thin curves assume A = 47r steradian. In this figure, we 
fixed Q = 16. 
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Figure 13. The galaxy power spectrum miiltiplicd by the mean number density, 
nPgai(fc, z = 0.9), as a function of A:. Here, wo fixed n{h/2f = 5x lO"-* [h-^M^c]-^ . We 
have nPgai(A:, z = 0.9) > 1 for < 0.2 /iMpc^^, where the shotnoise is subdominant. 



